First order quantum phase transitions 
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Quantum phase transitions have been the subject of intense investigations in the last two decades 
[1]. Among other problems, these phase transitions are relevant in the study of heavy fermion 
systems, high temperature superconductors and Bose-Einstein condensates. More recently there is 
increasing evidence that in many systems which are close to a quantum critical point (QCP) different 
phases are in competition. In this paper we show that the main effect of this competition is to give 
rise to inhomogeneous behavior associated with quantum first order transitions. These effects are 
described theoretically using an action that takes into account the competition between different 
order parameters. The method of the effective potential is used to calculate the quantum corrections 
to the classical functional. These corrections generally change the nature of the QCP and give rise 
to interesting effects even in the presence of non-critical fluctuations. An unexpected result is the 
appearance of an inhomogeneous phase with two values of the order parameter separated by a first 
order transition. Finally, we discuss the universal behavior of systems with a weak first order zero 
temperature transition in particular as the transition point is approached from finite temperatures. 
The thermodynamic behavior along this line is obtained and shown to present universal features. 

PACS numbers: 75.10.Jm ; 75.30.Kz ; 03.67.-a 
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INTRODUCTION 

Quantum phase transitions have been intensively stud- 
ied in the last two decades From a pure theoretical 
curiosity it became a field of intense experimental activ- 
ity in different areas of condensed matter physics. The 
basic concept in this field is that of a quantum criti- 
cal point (QCP). This is an unstable fixed point which 
separates a phase with long range order from a disor- 
dered phase at zero temperature 0, A fundamental 
distinction between this type of critical point and that 
associated with thermal phase transitions is the special 
role that time plays as an additional dimension. This is 
explicitly manifested in the quantum hyperscaling law, 
2 — a = v(d + z) which relates the dimension of the sys- 
tem d to the usual exponents a of the singular part of 
the free energy density and v the correlation length ex- 
ponent 0, Q . The new feature here is the appearance of 
the dynamic exponent z that arises from the time direc- 
tions. This appears in the suggestive form of an effective 
dimension d e f f = d + z which in fact controls the charac- 
ter of the quantum fluctuations and has most important 
consequences. For many problems of interest in the lab- 
oratory d e f f turns out to be larger or equal to the upper 
critical dimension d c of the problem and consequently all 
critical exponents are well known. In this case knowledge 
of the dynamic exponent is sufficient to characterize the 
universality class of the quantum phase transition. 

The scaling form of the free energy density close to a 
QCP is given by @, 
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where g measures the distance to the QCP (g = 0). This 
expression allows to obtain the dominant thermodynamic 
behavior of the system in the vicinity of the QCP. The 
field h is that conjugated to the order parameter and the 
exponents a, (3, 7, v are the usual critical exponents re- 
lated by standard scaling laws 3] . Only the hyperscaling 
relation is modified as discussed above. For d e ff > d c 
there may be dangerously irrelevant interactions which 
influence the critical behavior []| , in particular along the 
quantum critical trajectory g — 0, T — > and determine 
the critical line of finite temperature phase transitions 
I 

However, recent works Q are showing that, this is 
not all there is about quantum phase transitions with 
deff > d + z. The region close to a QCP seems to be 
a turbulent zone where many phases compete. The in- 
tensity of magnetic fluctuations near a magnetic QCP 
provides an additional mechanism for pair formation that 
favors the appearance of superconductivity 7] . In Kondo 
lattices local Kondo fluctuations may interfere with the 
long range magnetic correlations close to the magnetic 
QCP ||- In field induced Bose Einstein transitions in 
magnetic systems, soft elastic modes can couple to the 
spin-wave excitations with effects on the critical behavior 
|9( . Soft modes due to generic scaling invariance may also 
change the critical behavior [l(| and gauge fluctuations 
can affect charged systems in the vicinity of their quan- 
tum phase transitions [Tlj . Common to all these cases is 
that in the region of the phase diagram close to a QCP 
there is a competition of different types of fluctuations. 
Our aim here is to study the effects of this competition. 
In condensed matter physics, the most well known ex- 
ample of fluctuations of another field interfering with a 
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phase transition is that of the electromagnetic field on 
the thermal superconducting transition The quan- 
tum equivalent of this effect is known in quantum field 
theory as the Coleman- Weinberg mechanism 12]. 

COLEMAN- WEINBERG MECHANISM IN 
CONDENSED MATTER 



Finally, we point out that the coupling A has disap- 
peared from Eq. © for the effective potential due to 
dimensional transmutation (l2|. 

COMPETITION BETWEEN 
SUPERCONDUCTIVITY AND 
ANTIFERROMAGNETISM 



In the solid-state version of the Coleman -Weinberg 
mechanism |lz) . we consider a superconductor at T = 
represented by a complex field (ipi,<fi2) coupled to the 
electromagnetic field 0, El ■ 
the model is given by, 

1 



The Lagrangian density of 



L 



(2) 



We are using h = c = 1 units and the indices fi, v run 
from to d = 3. In Eq. J5J space and time are isotropic 
and consequently the dynamic critical exponent z = 1. 
For a neutral superfluid (q — 0) the system decouples 
from the electromagnetic field and has a continuous, zero 
temperature superfluid-insulator transition at m 2 = 
(see Fig. |SJ). 

The method of the effective potential Q yields the 
quantum corrections to the action given by the La- 
grangian density of Eq. ©. At T — in the one loop 
approximation, the effective potential close to the tran- 
sition (to « 0) is given by [j| 



1 



2,„2 



V eff=2 m V 
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(3) 



where (ip) is the minimum of the effective potential. We 
can show |l4l 1 1 ol] that the condition for such a minimum 
to exist is Xl -C £, where £ is the coherence length and 
Xl the London penetration depth as usually defined for 
Ginzburg-Landau models. In this case we find Q that 
at a critical value of the mass, given by 



32tt : 



(4) 



there is a first order transition to a superconducting 
state. Notice that the transition in the neutral super- 
fluid (q = 0) is continuous rather than first order and 
takes place at to 2 = 0. Therefore, the coupling to the 
electromagnetic field in the charged superfluid lead to 
symmetry breaking, shifting the transition of the neu- 
tral superfluid (see Fig. |SJ) and changed its nature from 
continuous to first order. The shift of the transition Eq. 
(@J depends on the coupling of the order parameter to 
the soft modes, in the present case, the charge q of the 
Cooper pairs. 



In this section we show that weak first order quantum 
phase transitions (WFOQPT) and spontaneous symme- 
try breaking can also occur due to the competition be- 
tween different types of instabilities in the same region 
of the phase diagram. 

We consider a Ginzburg-Landau model which is appro- 
priate to describe the competition between superconduc- 
tivity (SC) and antifcrromagnetism in a heavy fermion 
metal. The model contains three real fields. Two fields, 
4>\ and fa, correspond to the two components of the su- 
perconductor order parameter. The other field fa, for 
simplicity represents a one component antifcrromagnetic 
(AF) order parameter. The free functional of the mag- 
netic part [16] takes into account the dissipative nature 
of the paramagnons near the magnetic phase transition 
in the metal [2j and is associated with the propagator, 



A)(w,q) 



i\uj\t - q 2 - to 2 



(5) 



where r is a characteristic relaxation time and to 2 gives 
the distance to the magnetic transition. The quadratic 
form of the superconductor, the same used in the previ- 
ous section, is given by, 



Go(fc) = G (^,q) = 



k 2 — to 2 



(6) 



The part of the action associated with the classical po- 
tential is, 

+V s (ct>i, fa) + V p (fa) + Vi(fa,fa,fa), (7) 
where the self-interaction of the superconductor field is 

A 



and that of the antiferromagnet, 

v P {fa) = |4 



(8) 



(9) 



Finally, the last term is the (minimum) interaction be- 
tween the relevant fields, 



Vi{fa,fa,fa) = u{<t>t + 44)44 



(10) 



This term is the first allowed by symmetry on a series ex- 
pansion of the interaction. Notice that for u > 0, which is 
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the case here, superconductivity and antiferromagnetism 
are in competition and this term does not break any sym- 
metry of the original model. However, including quan- 
tum fluctuations we show that spontaneous symmetry 
breaking can occur in the normal phase separating the 
SC and AF phases. 

The first quantum correction to the potential can 
be obtained by the summation of all one loop dia- 
grams (Fig.HJ. 
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FIG. 1: 

tential. 



One loop diagrams contributing to the effective po- 



We apply the general method proposed by Cole- 
man |l2l | with minimum modifications to account for the 
different nature of the propagators in our problem. The 
sum over the field indices can be easily done if we define 
a vertex matrix M, given by 

[«]„„ - 



(11) 

and then take the trace. In Eq. (f 1 1 1> the propagators 
(.Kg = Go or Do) are incorporated in the definition of the 
matrix. We draw the loops with arrows and choose the 
outgoing propagator of each vertex to be included in the 
associated element. The matrix M is then obtained de- 
riving the classical potential with respect to the fields {(/>} 
and taking the values of these derivatives at the classical 
values of the fields, {4>i C }- The sum of diagrams with 
the correct Wick factors is formally done in momentum 
space and using the property of the trace 



Tr[ln(l - M)} = lndet[(I - M)], 



we get 



V w [(t> c ] = l -h J d 4 fclndet [I 



M(k)] 



(12) 



(13) 



The 3x3 matrix M can be simplified if we choose the 
classical minimum of the superconductor fields imposing 
4>2c = (this can be done because the minimum depends 
only on the modulus 4>\ c + 0|c)- Hence, rotating to Eu- 
clidean space, so that, k 2 = uj 2 + q 2 and using h = 1 units 
the first quantum correction can be written as 



V (1) (01c> 



In 
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k 2 + m 

C(01c,03c) 



>lc, <P3c 



(k 2 + m 2 )(\u)\T + q 2 + m 2 ) 



(14) 



A^ lc ,cf> 3c ) = (X/6)cf> 2 lc + 2u4 c 



B(<t>i c ,4> 3c ) = (A/2)0? C + 2 U 



3c 



C(& CJ &c)=2u0? c +(ff/2)$ o 



-D(01e,03c) = 4u0i c 3c 



(15) 
(16) 
(17) 
(18) 



The total effective potential with first order quantum cor- 
rections is then given by 



eff\ 



he, 03c) = Wlc, 03c) + ^ (1) (01c, 03c) (19) 



where V c i is the classical potential of Eq. 10 and V^ 1 ' is 
the first quantum correction of order h of Eq. (|14(l . 



Superconducting transition 

We first consider the effect on the superconductor tran- 
sition in HF in the presence of antiferromagnetic para- 
magnons (0i c ^ 0, 03 C = 0). Detailed calculation of the 
effective potential have already been presented 0] . The 
general result is given by, 



Veff(<l>. 



+am>Ll0ic| + ^0ic+O(0 5 )- (20) 



where 



In Eq. H2Ufl . M is a renormalized superconducting mass 
and A a renormalized coupling, of the same order of the 
bare coupling A. The new coupling a introduced by fluc- 
tuations can produce a symmetry breaking in the normal 
state extending the SC region in the phase diagram at 
T = 0. The same mechanism turns the superconduct- 
ing transition, which was continuous before coupling to 
the paramagnons, to weak first order with a small latent 
heat 0| . Introducing again the coherence length £ and 
the London penetration depth A/, we can show that the 
condition for the existence of minima away from the ori- 
gin is equivalent to Al £ as in the previous case [lij . 
It's also interesting to notice that the first order transi- 
tion is produced by the cubic term in Eq. l|2UI) and this 
term is proportional to the magnetic mass m p . If mag- 
netic fluctuations were critical, i.e., m p = 0, the only ef- 
fect of the coupling would appear as a term proportional 
to 5 . This term is usually neglected since its power 
is higher than those initially considered in the classical 
potential and usually insufficient to create new minima 
around the origin. Therefore, if AF fluctuations were 
critical the effects of the quantum corrections in the tran- 
sition could be neglected. 
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Antiferromagnetic transition 

We now study the effect of SC fluctuations in the mag- 
netic transition (<f> lc — 0, foe ^ 0). We consider two 
kinds of quadratic forms associated with the free super- 
conducting fields. The first is the usual Lorentz invariant 
free action used before. Next, we work with another free 
action which takes into account dissi pat ion and is asso- 
ciated with a z = 2 dynamics [ItI fisl Il9l| . 

Let's first study the case of the Lorentz invariant prop- 
agator. Close to the AF transition we have 



Veffifoc) « -MX + |<4 + u 2 4c In 



+um 2 (f)l In 



(21) 



Notice from Eq. I|21|) that if the SC fluctuations were 



critical, i.e. 



= 0, wc would obtain the same result 



of Eq. Q, i.e., a fluctuation induced quantum first or- 
der transition. However, if the SC fluctuations are not 
critical but close to criticality the last term of Eq. i|21|) 
may become important. Of course, its relevance depends 
on the strength of the renormalized coupling u and the 
results considering this term lead to new and interest- 
ing changes in the ground state. We obtain, besides the 
two finite minima of the Coleman- Weinberg potential of 
Eq. @, two extra minima very close to the origin |l6| . 
The states associated with these minima have a small 
value of the order parameter, the sub-lattice magnetiza- 
tion. An additional first order transition occurs when 
the other two minima away from the origin become the 
stable ones as the system moves away from the supercon- 
ductor instability. This transition is from a small moment 
AF (SMAF) to a large moment AF (LMAF) and occurs 
even before the continuous mean field transition. When 
we move away from the magnetic transition, i.e., towards 
the superconducting instability, the strength of this new 
term decreases with the value of m 2 and the two new min- 
ima move to the origin producing a normal state with 
vanishing sub-lattice magnetization again. Notice that 
the SMAF phase is obtained because the magnetic order 
parameter couples to superconducting fluctuations which 
are non-critical. Critical fluctuations yield the same re- 
sults of section B 

Now, for many cases of interest, SC fluctuations are 
better described by a dissipative propagator associated 
with az=2 dynamics 0,^3 similar to Eq. JSJ. This is 
useful to account for pair breaking interactions, as mag 
netic impurities that can destroy superconductivity 
It is given by, 



G (w,q) 



i|w|r' — q 2 — m 2 



(22) 



The parameter m 2 is still related to the distance from 
the SC phase transition and we have a relaxation time 




V 



cf 




FIG. 2: New minima appear in the potential for urn 2 ^ 0. 
The effective potential is shown here for two situations: for 



Mp , where the first order transition from LMAF to 



SMAF occurs and these two states become degenerate and 
for the spinodal point at which the LMAF becomes unstable 
inside the SMAF phase. 



t' . Calculation of the effective potential is very similar 
to the previous cases and the result has the form 



Veff = ~ 2 M 2 



4!' 



+ 



— (2u4>i+m 2 fl 2 (23) 



with a renormalized magnetic mass M p and coupling g. 
Quantum corrections can once again produce a weak first 
order transition. An analysis of the extrema of the effec- 
tive potential Eq. (|23[l shows that the transition can be 
first order depending on the coupling values. The ap- 
pearance of SMAF phases is not possible in this case. 



COUPLING TO LOCAL MODES 

We now study the coupling of antiferromagnetic (AF) 
fluctuations to local modes in a three dimensional sys- 
tem. This model is useful to describe the QCP of heavy 
fermions where local modes can coexist with antiferro- 
magnetic fluctuations. The local modes can be either 
Kondo Q or valence fluctuations [53] • The local propa- 
gator in Euclidean space is written as 



1 



M/7 



(24) 



where m 2 gives the distance to the QCP and r is asso- 



ciated with the lifetime of the excitations, 
paramagnons, we have 



G v {u,q)=G p (k) 



1 



For the AF 



(25) 
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The action is 

s = J d^G-'ikM 2 



12 



M 4 + g-» 



(26) 



Let us consider one component fields, the classical poten- 
tial (not including the mass terms) is given by, 



(27) 



and the effective potential by, 
1 f d 4 fc 



In 1 



2 J (2tt) 4 
d 4 k (. gift 



(27T) 



1 



H 



(28) 



where the dots represent counter-terms |21|. The inte- 
gration in u> in the first integral can be performed us- 
ing a cut-off, A z — m 2 L T (z = 2). However, the integral 
over <iq diverges when A — * oo and this term is non- 
renormalizable unless we introduce a cut-off q c and for 
later purposes f = (Aq^)/(3n). 



Case m 2 L 3> fu 2 



In this case the analysis of the renormalized effective 
potential shows that, as for the classical potential, there 
is a second order phase transition at to 2 = and the 
coupling to the local modes has no significant effects. 



Case m\ < fu 2 

This is the most interesting case. The renormalized ef- 
fective potential after dimensional transmutation is given 
by, 



-fuip In 



V> 2 



8 / 6m 2 7r 2 + fu 



5/2 
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(V>) 2 7T 2 



ij) 2 + to 2 — fm 2 L ln(uip 2 ) 



(29) 



where the renormalized quantities M 2 and G are inde- 
pendent of g [5lJ ■ In this case we find there is no tran- 
sition for m 2 > 0, such that, there is no dynamic sym- 
metry breaking. A first order transition occurs for neg- 
ative values of m 2 , at a critical value vn c2 which has to 
be obtained numerically. The limit of stability of the 
minima of the effective potential (spinodal) is given by, 

1 For the par- 



m'p Sp = — (fu)/(47r 2 ), with m 
ticular case to 2 



2 > m c< 

p,sp ^ p 

there is a first order transition be- 



tween the paramagnetic and antiferromagnetic phases, as 



shown in Fig. [3J Then, even in 3d, when both transitions 
coincide, the magnetic transition becomes first order due 
to the effect of the critical local fluctuations. Now, for 




FIG. 3: For m L = there is a first order transition to the 
antiferromagnetic phase [2l|. 

to 2 small but ^ 0, the phase diagram changes drastically 
since new minima appear in the potential close to the ori- 
gin. These minima which correspond to small values of 
the order parameter (SMAF) can coexist with those as- 
sociated with the large values of this parameter (LMAF). 
As m 2 further decreases there is a first order transition 
between the SMAF and LMAF phases (see fig. EJ). 



3.0x1 0^ 

Vfo) 




FIG. 4: 
gin [p. 



For m 2 L 7^ new minima appear close to the ori- 



SCALING AT A WEAK FIRST ORDER 
QUANTUM TRANSITION 

At a first order transition there is no true critical 
behavior since the correlation length does not diverge. 
However it turned out to be useful to develop a sc aling 
approach for these transitions in the classical case [22j. 
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As we show below the same is true for first order quantum 
phase transitions. This is particularly useful for WFO- 
QPT where we expect the correlation length and char- 
acteristic time to become very large. The best way to 
introduce these ideas is to consider a specific case, for 
example, the Coleman- Weinberg transition of section D 

Introducing a parameter g = to 2 — m 2 which measures 
the distance to the first order transition at m c , we can 
write Eq. |2| for the effective potential at T — and for g 
small as, 



eff 



to 2 — to 2 | oc \g\ 2 ~ 



(30) 



with the exponent a = 1 reflecting the fact the tran- 
sition is first order |l3l| . The associated latent heat is 
Lh = \m 2 U£^. Spinodal points at T = can also be 
calculated 1131. 



T 



Line of First Order , 
Transitions 
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the shift exponent ip" 1 — z/(d + z — 2) = 1/2 in d = 3 
(see Ref. |3j). The new line of first order transitions is 
given by T c oc \Jm 2 — |m 2 |. 

We will now consider the system sitting at the new 
quantum phase transition point, i.e., at m = to 2 and de- 
crease the temperature. For high temperatures, T 3> m c , 
which corresponds to the regime I of Fig- EJ the function 
hiv) saturates, I%{y < 0.12) w —2.16. In this case the 
effective potential, 



V eff {T)^\m 2 ^) 2 \g\\\ 



4.32 T d+1 



ir 2 m 2 (<p) 2 
and can be cast in the scaling form, 



1.91 



Veff(T)<x\g\ 



"F 



T 



with F(0) = constant. This scaling form is reminiscent 
of that for the free energy close to a quantum critical 
point. In the present case of a discontinuous zero tem- 
perature transition, the critical exponent 13] a = 1 and 
the characteristic temperature is, 



« |<?r ex \g\ 



1.91 



1.91 



with^ = l/(d+z) [H. In this regime I or scaling regime, 

along the line m 2 = m 2 shown in Fig. the free energy 

density has therefore the scaling form f(m — m c ,T) oc 
rp(d+z)/z anc i ^ e S p ec ig c hgg^ j s gi ven by, 



m 



C/T\, „, oc T 

I \(m—m c ,Tj 



(32) 



FIG. 5: Phase diagram of a charged superfluid coupled to 
photons. For completeness we show also the critical line of 
the neutral superfluid. Along the trajectory m 2 = m 2 one 
can distinguish different regimes as explained in the text. 



The finite temperature case can be studied replacing 
the frequency integrations in the calculation of the ef- 
fective potential by sums over Matsubara frequencies Q . 
The effective potential at finite temperatures close to the 
transition can be written as Il3j 



V e}f (T) = \ m 2 { V ) 2 \g\ 



1- 



TT 2 m 2 (ip} 2 
where the integral / is given by 



2 T d+1 fM(cp) 



id(y) 



dxx d 1 In 



1 



(31) 



and M(ip) = m 2 + q 2 tp 2 . The function ^(y) = I(y) can 
be obtained numerically integrating Eq. I|31|l . 

The finite temperature phase diagram is shown in 
Fig. 03 For completeness we show the critical line of the 
neutral superfluid, T$f oc |m 2 |^, which is governed by 



Then the thermodynamic behavior along the line m 2 = 
m 2 in regime I (T> m c ) is the same as when approach- 
ing the quantum critical point of the neutral superfluid, 
along the critical trajectory m 2 = 0. The system is un- 
aware of the change in the nature of the zero temperature 
transition and at such high temperatures charge is irrel- 
evant. 

When further decreasing temperature along the line 
m 2 = m 2 there is an intermediate, non-universal regime 
(regime II in Figs.|SJ). In the present case for T ~ m c the 
specific heat C/T d ' z oc InT 13]. 

Finally, at very low temperatures, for T « m c and 
m 2 = rn 2 , i.e., in regime III of Fig. the specific 
heat vanishes exponentially with temperature, C ' jT d l z oc 
exp(— m c /T). The gap for thermal excitations is given by 
the shift m c of the quantum phase transition. The cor- 
relation length which grows along the line to 2 = to 2 with 
decreasing temperature reaches saturation in regime III 
at a value £s which is function of the inverse of the gap. 
The exponential dependence of the specific heat is due 
to gapped excitations inside superconducting bubbles of 
finite size £s oc m c ■ 

Although the results above have been obtained for the 
model of section D the behavior in the scaling regime I 
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and III should be universal and characteristic of any weak 
first order quantum transition. Notice that the relevant 
critical exponents which determine the scaling behavior 
in particular in regime I are those associated with the 
QCP of the uncoupled system which in the present case 
is the neutral supcrfluid. 
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